Scaling the localisation lengths for two 
interacting particles in one-dimensional 
random potentials 

Rudolf A. R6mer a , Mark Leadbeater b , and Michael Schreiber a 

a Institut fur Physik, Technische Universitdt, D-09 107 Chemnitz, Germany 

b Dipartimento di Fisica, Universita di Roma Tre, Via della Vasca Navale 84, 

I- 00 146 Roma, Italy 

Abstract 

Using a numerical decimation method, we compute the localisation length A2 for 
two onsite interacting particles (TIP) in a one-dimensional random potential. We 
show that an interaction U > does lead to \2(U) > A2(0) for not too large U and 
test the validity of various proposed fit functions for A2(i7)- Finite-size scaling allows 
us to obtain infinite sample size estimates ^(U) and we find that &(U) ~ £2(0) 
with a(U) varying between a(0) 1 and a(l) 1.5. We observe that all ^(U) 
data can be made to coalesce onto a single scaling curve. We also present results for 
the problem of TIP in two different random potentials corresponding to interacting 
electron-hole pairs. 



In two recent articles [1,2], we studied as a simple and tractable approach to 
the problem of interacting electrons in disordered materials the case of only two 
interacting particles (TIP) in ID random potentials. Previous considerations 
[3] had led to the idea that attractive as well as repulsive interactions between 
TIP give rise to the formation of particle pairs whose localisation length A2 is 
much larger than the single-particle (SP) localisation length Ai ~ 105 /W 2 , 

A 2 ~ U 2 \\ (1) 



at two-particle energy E — 0, with U the Hubbard interaction strength. Al- 
though many papers have numerically investigated the TIP effect [3-9], an 
unambiguous reproduction of Eq. (1) is still lacking. However, it appears well 
established that some TIP delocalisation such as A2 > Ai does indeed exist 
due to the interaction. Recently, a duality in the spectral statistics for U and 
y/24/U has been proposed [11] for small and very large \U\. 
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Fig. 1. X 2 (U) for TIP as a function of \U\ (solid lines) and \/24/|J7| (dashed lines) 
at E = for disorders W = 3 (+), 4 (*), and 5 (x) and M = 201. The data are 
averaged over 100 samples. The lines (symbols) indicate data for U > (U < 0). 

In Refs. [1,2], we have employed a numerical decimation method [10], i.e., 
we replaced the full Hamiltonian by an effective Hamiltonian for the doubly- 
occupied sites only. In [1], we considered the case of TIP with n, m corre- 
sponding to the positions of each particle on a chain of length M and random 
potentials e l n = e 2 n G [— W/2, W/2]. In [2], we studied the case where e l n and 
are chosen independently from the interval [— W/2, W/2], which may be 
viewed as corresponding to an electron and a hole on the same chain (IEH). 
Via a simple inversion, we then obtained the Green function matrix elements 
(1, 1\G 2 \M, M) between doubly-occupied sites (1, 1) and (M,M) and focused 
on the localisation length A2 obtained from the decay of the transmission 
probability from one end of the system to the other, i.e., 

i = _ ^~ii ln|(1 ' 1|G2|M ' M)| - (2) 

In Fig. 1 we present data for \2(U) obtained for three different disorders for 
system sizes M = 201 at E = 0. In agreement with the previous arguments 
and calculations [6,7,11], we find that the enhancement is symmetric in U and 
decreases for large \U\. In [11] is has been argued that at least for Ai ~ M, 
there exists a critical U c = 24 1 / 4 w 2.21, which should be independent of W, 
at which the enhancement is maximal. We find that in the present case with 
Ai < M the maximum of A 2 (£/) depends somewhat on the specific value of 
disorder used. The data in Fig. 1 may be compatible with the duality of Ref. 
[11], but only for the large disorder W = 5. For the smaller disorders and for 
the range of interactions shown, we do not observe the duality. We emphasize 
that the duality observed in [11] is for spectral statistics and need not apply 
to quantities such as the localisation length A2. 
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In order to reduce the possible influence of the finiteness of the chain length, 
we constructed finite-size-scaling (FSS) curves for 11 interaction values U = 
0, 0.1, . . . , 1 from the A 2 data for 26 disorder values W between 0.5 and 9, for 24 
system sizes M between 51 and 251, averaging over 1000 samples in each case. 
In Fig. 2 we show the infinite-size localisation lengths (scaling parameters) £ 2 
obtained from these 11 FSS curves. A simple power- law fit £ 2 oc W~ 2a in the 




Fig. 2. Left panel: TIP localisation lengths £2 after FSS. The dashed lines repre- 
sent power-law fits. Inset: Exponent a obtained by the power-law fits. Right panel: 
Scaling plot according to [5] with TIP localisation lengths &(U) for W 6 [1,5]. The 
solid line indicates a slope of 1/4, the dashed line the value of in the limit 

U = 0. 

disorder range W G [1,5] yields an exponent a which increases with increasing 
U as shown in the inset of Fig. 2, e.g., a = 1.55 for U — 1 and a = 1.1 for 
U — 0. Because of the latter, in the following we will compare ^ 2 (U 7^ 0) with 
£2(0) when trying to identify an enhancement of the localisation lengths due 
to interaction. 

Song and Kim [5] suggested that the TIP localisation data may be described 
by a scaling form £ 2 = W~ a °g(\U\/W A ) with g a scaling function. They obtain 
A = 4 by fitting the data. Our data can be best described when a is related 
to the disorder dependence of £2 as (a — eto)/A ~ 1/4. As shown in Fig. 2, 
the scaling is only good for W G [1,5] and U > 0.3. We note that assuming 
an interaction dependent exponent a(U), we still do not obtain a good fit to 
the scaling function with the data for all U. 

In Fig. 3, we show that a much better scaling can be obtained when plotting 

UU)-U0) = g[f(UM0)] (3) 



with f(U) determined by FSS. Now the scaling is valid for all U and W G 
[0.6, 9]. As indicated by the straight lines, we observe a crossover from a slope 
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2 to a slope 3/2. There are some deviations from scaling, but these occur 
for large and very small values of 6(^0 and are most likely due to numerical 
inaccuracy [1]. In the inset of Fig. 3, we show the behavior of f{U). For U > 0.3 
a linear behavior f(U) oc U appears to be valid which translates into a U 2 
(U 3/2 ) dependence of &(U) - 6(0) in the regions of Fig. 3 with slope 2 (3/2). 
For U < 0.5, we have f(U) oc VU which yields &(U) - 6(0) oc U (f/ 3/4 ). 




f(u)%(0) m%(0) 



Fig. 3. Left panel: Scaling plot of Eq. (3) for TIP with W G [0.6,9]. The solid 
(broken) line indicates the slope 2 (1.5). Inset: The values of f(U) needed to make 
the data collapse onto the U = 0.1 curve. Right panel: IEH localisation length scaled 
as in the left panel. The solid line indicates slope 1.61 and W £ [1, 7]. 

Thus in summary it appears that our data cannot be described by a simple 
power-law behavior with a single exponent as in Eq. (1) neither as function of 
W, nor as function of 6(0) nor after scaling the data onto a single scaling 
curve. 

As for TIP we computed [2] the IEH localisation lengths by the DM along the 
diagonal using 100 realizations for each (U, M, W). We find that the data for 
IEH are very similar to the case of TIP. We again perform FSS and observe that 
the infinite-size estimates 6(^0 are well characterized by an exponent a(U). 
We can again scale the £(U) data for IEH onto a single curve as shown in Fig. 
3. However, here the crossover from slope 2 to 3/2 is much less prominent and 
the data can be described reasonably well by a single slope of 1.61. Also, the 
crossover behavior in f(U) is suppressed. We remark that these differences 
may be due to the smaller number of samples used for IEH. 

In conclusion, we observe an enhancement of the two-particle localisation 
length due to onsite interaction both for TIP and IEH. This enhancement 
persists, unlike for TMM [6,8,9], in the limit of large system size and after 
constructing infinite-sample-size estimates from the FSS curves. We remark 
that the IEH case is of relevance for a proposed experimental test of the TIP 
effect [12]. 
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